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INTRODUCTION 

Robust  prediction  of  the  motion  of  debris 
from  explosions  is  critical  to  the  estimation  of 
separation  distances,  fragment  density,  and 
debris  lethality  in  explosive  safety  hazard 
analysis.  Debris  dispersion  models  have  been 
developed  and  applied  to  many  specific  types 
of  hazardous  debris  (ranging  from  masonry, 
concrete,  structural  steel  elements,  soil  and 
rock  ejecta,  to  bomb  fragments).  These 
models  generally  use  2-D  (degree-of- 
freedom)  or  3-DOF  trajectory  calculations 
and  consider  only  the  drag  component  of  the 
aerodynamic  force  vector  (e.g.,  DOE/TIC- 
11268  [1980],  Huang  [1984],  McCleskey 
[1988],  Bowles  and  Oswald  [1990]).  For 
many  debris  geometries  (such  as  plates, 
slender  fragments,  and  structural  elements 
with  high  slenderness  ratios),  drag  component 
trajectory  models  underpredict  maximum 
debris  range  and  dispersion.  In  addition,  if 
the  potential  for  spin- stabilized  motion  exists, 
the  safety  distance  for  plate  or  disk  shape 
fragments  can  significantly  exceed  that  which 
would  be  predicted  for  random  tumbling  or 
other  non-spin-stabilized  motion  [Twisdale, 
1984], 

This  paper  presents  several  topics  related 
to  explosive  safety  hazardous  debris  trajectory 
analysis.  The  scope  of  the  paper  is  limited  to 
free-flight  trajectory  analysis  and  does  not 
include  debris  ricochet  or  ground  roll  models. 
Following  a  review  of  alternative  trajectory 
models,  key  features  of  the  random 
orientation  (RO)  6-D  model  are  summarized. 


The  RO  6-D  model  uses  drag,  lift,  and  side 
force  components  and  randomly  updates  the 
rigid  body  orientation  of  the  missile.  It  has 
been  implemented  in  the  TORMIS  and 
TURMIS  computer  codes  [Twisdale,  et  al., 
1978,  1979,  1981,  1984]  for  facility  risk 
assessment  from  wind-borne  debris  hazards 
and  for  fragments  and  secondary  missiles 
from  exploding  equipment.  An  aerodynamic 
library  from  existing  databases  has  been 
developed  for  fourteen  generic  missile/debris 
shapes  (including  structural  components, 
plates,  chunky  fragments,  and  secondary 
missiles  from  internal  equipment).  An 
analytically  derived  random  tumbling  mode 
drag  coefficient  (as  a  function  of  the  axial  and 
cross-flow  coefficients)  for  cylindrical  shapes 
is  summarized  for  use  with  2-D  transport 
predictions  of  chunky  fragments.  An 
equation  is  also  presented  that  allows  an 
evaluation  of  whether  or  not  spin-stabilized 
motion  will  occur  for  in-plane  rotation  of 
discs  and  plates.  This  paper  concludes  with 
RO  6-D  vs.  2-D  drag  comparisons  of  free- 
flight  debris  range  prediction  for  several 
secondary  debris  shapes. 

REVIEW  OF  TRAJECTORY  MODELS 

Debris  transport  methodology  predicts  the 
free-flight  motion  of  the  primary  fragments 
and  secondary  missiles  that  are  generated  by 
the  explosion.  A  set  of  initial  conditions  are 
required  for  the  trajectory  analysis,  including: 
the  missile  debris  mass,  geometry,  initial 
translational  and  angular  velocities,  ejection 
angles,  and  missile  inertial  orientations. 
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Given  these  initial  conditions,  the  transport 
methodology  consists  of  aerodynamic  models 
of  the  missile  shapes,  the  governing  dynamic 
and  kinematic  relations,  and  the  solution 
scheme  for  the  developed  equations  of 
motion.  Integration  of  these  equations  yields 
the  motion  time- history  of  the  missile,  which 
provides: the  means  to  predict  the  free-flight 
motion,  impact  conditions  and  density, 
lethality,  and  safety  distance. 

Table  1  summarizes  several  basic 
transport  models  that  are  available  for 
explosive  safety  debris  hazard  analysis. 
Trajectory  models  are  most  commonly 
distinguished  by  the  type  of  motion  they 
describe.  Generally,  one  degree  of  freedom 
(1-D)  refers  to  motion  of  a  point  along  a  line; 
two  degree  of  freedom  (2-D),  to  a  point  in  a 
plane;  three  degree  of  freedom  (3-D),  to 
motion  of  a  point  in  space;  and  six  degree  of 
freedom  (6-D),  to  translational  and  rotational 
motion  of  a  rigid  body  in  space.  Another 
distinguishing  feature  is  the  number  of 


aerodynamic  force  components,  including 
moments,  that  are  considered. 

The  simplest  model  in  Table  1  is  the  2-D 
(2  degree  of  freedom)  model  for  a  particle 
mass  subjected  only  to  the  force  of  gravity. 
Two  ordinary  differential  equations,  which 
can  be  integrated  in  closed  form,  describe  the 
parabolic  motion  of  the  particle  within  a 
vertical  plane.  This  model  generally  is  valid 
only  forshort  distance  trajectories  within  the 
donor  facility  to  get  impact  conditions  for 
sympathetic  detonations  and/or  secondary 
debris  generation.  The  next  hierarchy  of 
model  sophistication  involves  the  introduction 
of  an  aerodynamic  drag  force  in  the  ballistic 
2-D  model.  The  advantage  of  this  model  over 
the  no  drag  model  is  that  it  provides  much 
more  accurate  predictions  of  motion,  impact 
speed,  and  position.  This  model  is  valid  for 
non-spinning  spherical  and  chunky  debris,  for 
which  lift  and  side  forces  and  moments  are 
negligible.  The  resulting  coupled  ordinary 
differential  equations  are  integrated 


TABLE  1.  SUMMARY  OF  ALTERNATIVE  TRANSPORT  MODELS. 


Features 

2-D,  No  Drag 

2-D.  Drag 

Models 

3-D 

RO  6-D 

6-D 

Parameters 

8 

8 

8 

8.  v 

/ 

Aerodynamic  Forces 

None 

cD 

-  Go 

cD,  CL ,  cs 

cD ,  CL,  CS, 

Cl  Cm,  Cp 

Equations  of  Motion3 

2  ODE 

2  Coupled 
ODE 

3  Coupled 
ODE 

3  Coupled 
ODE 

3  Force  Eq. 

6  Coupled 
ODE 

SimulationEfficiency 

Analytic 

High 

High 

Moderate 

Low 

Impact  Speed  Prediction*3 

+ 

~ 

~ 

Yes 

Yes 

Impact  Position  Predictionb 

+ 

- 

- 

- 

Yes 

Impact  Dispersion 

- 

- 

- 

Yes 

Impact  Orientation  Predictionb 

No 

No 

No 

~ 

Yes 

Impact  Obliquity  Prediction*3 

~ 

~ 

~ 

Yes 

Yes 

Impact  Angular  Velocity  Prediction*3 

No 

No 

No 

No 

Yes 

a  ODE  =dDrdiriary  Differential  Equations. 

b  ~  =  Approximately  Correct 
-  =  Tendency  to  Underestimate 
+  =  Tendency  to  Overestimate 
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numerically  to  predict  debris  motion  time 
history.  The  3-D  model  in  Table  1  predicts 
the  general  motion  of  a  point  in  space.  Its 
basic  parameter  is  also  the  drag  coefficient, 
whose  value  is  often  specified  to  account  for 
random  tumbling  of  the  object.  The  6-D 
models  in  Table  1  simulate  the  aerodynamics 
of  rigid  bodies  that  cannot  be  adequately 
treated  by  the  simpler  2-D  and  3-D  models. 
The  random  orientation  model  (RO  6-D) 
considers  drag,  lift,  and  side  forces  and 
simulates  missile  tumbling  by  periodic 
reorientation  [Twisdale,  1979],  Its  prediction 
capabilities  are  enhanced  over  the  particle 
models  with  only  a  modest  decrease  in 
simulation  efficiency.  Conventional  6-D 
models  [Etkin,  1977;  Redmann,  et  al.,  1978] 
track  missile  translation  and  rotation  using  a 
system  of  six  coupled,  ordinary,  nonlinear 
differential  equations.  Such  models  require 
estimation  of  aerodynamic  force  and  moment 
coefficients  over  all  body  orientations,  which 
are  generally  not  known  for  arbitrary  bluff- 
body  fragments. 

RANDOM  ORIENTATION  MODEL 

This  section  presents  the  key  features  of 
the  RO  6-D  model  for  debris  dispersion 
analysis,  summarized  from  Twisdale,  et  al. 
[1978,  1979,  1981,  1984],  In  this  model,  the 
actual  rigid  body  orientation  of  the  missile  is 
considered  and  the  aerodynamic  specification 
includes  drag,  lift,  and  side  force  components. 
Figure  la  shows  the  inertial  reference  frame 
along  with  other  reference  frames  that  will  be 
used  in  the  development  of  the  model. 

The  missile  centerline  orientation  is 
specified  by  two  randomly  determined  angles 

(y/,  (p )  measured  from  a  (u,  v,  w)  coordinate 
system  as  defined  in  Figure  lb.  The  relative 

velocity  vector  defines  the  v  direction,  while 
u  =  (O  x  k)/\  0  x  k\  and  w  =  {u  x  0).  Once  the 
missile  orientation  is  established  for  a  time 
step,  wind  axis  unit  vectors  are  determined  by 
forming  the  vector  cross  product  of  the 
missile  centerline  position  unit  vector  ( MCL ) 


a.  Inertial  Reference  Frame 


b.  Missile  Orientation 


c.  Relative  Wind  Frame 


Figure  I,  Coordinate  Systems  and  Missile 
Specification. 
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with  the  relative  velocity  vector  ( v )  to 
establish  the  pitch  axis  ( P ).  The  missile 
diameter  unit  vector  (MD)  is  rotated  through 
a  randomly  selected  angle  (5)  from  the  pitch 
axis.  The  relative  velocity  unit  vector  fu)  is 
then  combined  with  the  pitch  axis  ( P )  in  a 
vector  cross  product  to  establish  the  lift  unit 
vector  (L).  This  approach  defines  the  wind 
axis  System  ( v,  P,  and  L)  for  each  time  step 
and  provides  the  respective  directions  for  the 
three  aerodynamic  force  components  (drag, 
lift,  and  side)  properly  oriented  for  the  missile 
attitude. 


The  magnitudes  of  the  three  translational 
forces~~sre  taken  as  proportional  to  the  square 
of  the  relative  velocity  and  to  the  three 
aerodynamic  coefficients  (Co,  C$,  and  C/J, 
which  may  each  be  functions  of  total  wind 

angle  of  attack  (a)  and  roll  angle  (5).  These 
angles  (a,  8)  are  both  shown  in  Figure  1c. 
The  missile  angles,  and  hence  vectors  MCL 
and  MD,  are  updated  at  selected  intervals 
according  to 


i//=ca$~l  (1  -2^) 

<P  =  K{ 2£2 -1) 

cx  —  cos-'  (sin  i/rcos  <p) 

S  -  2  7T^2 


0  <  y/<  k 
-n<(j)  <jc 
0  <  a  <  n  _ 
0  <  8  <  2  n 


(1) 


where  §j,  ^2.  and  £3  are  random  numbers 
selected  from  a  uniform  distribution  on  the 
unit  interval.  The  time  between  missile 
orientation  updates  is  termed  the  update 
period,  and  its  reciprocal,  update  frequency. 

The  angles  a  and  8  are  used  as  input  to  the 
aerodynamic  coefficient  determination.  Once 
the  three  coefficients  are  determined,  they  are 
combined  with  the  dynamic  pressure, 
reference  area  (/l),  and  the  three  appropriate 
wind  axis  unit  vectors  to  form  the  total 
aerodynamic  force  for  a  single  time  step. 

The  mass  center  of  the  missile  is  tracked 
relative  to  the  reference  frame,  according  to 


the  standard  dynamic  equations  of  motion. 
These  equations  form  a  set  of  six  coupled, 
nonlinear,  ordinary  differential  equations  that 
define  an  initial  value  problem  for  a  set  of 
prescribed  initial  conditions.  Shampine’s 
method  [Shampine  and  Gordon,  1975]  is  used 
to  integrate  these  equations. 

MISSILE  AERODYNAMICS  FOR 

RO  6-D  MODEL 

Since  complete  aerodynamic 
characteristics  generally  do  not  exist  for  the 
potentially  wide  variety  of  debris  shapes,  a 
modified  cross-flow  theory  has  been  applied 
to  develop  the  aerodynamic  coefficients  for 
the  random  orientation  model.  This  approach 
has  been  successfully  used  to  develop  the 
wind  axis  aerodynamic  forces  as  a  function  of 
angle  of  attack  for  slender  cylinders  knowing 
only  the  drag  force  coefficients  for  the  body 
in  normal  flow  to  the  major  body  axes 
[Hoemer,  1965]. 

The  basic  theory  assumes  the 
superposition  of  two  flows  perpendicular  to 
the  missile  axis  (axial  and  cross  flow)  in 
which  the  magnitude  of  the  mutually 
orthogonal  flows  is  determined  vectorally 
knowing  freestream  velocity  and  angle  of 
attack.  The  aerodynamic  forces  acting  on  the 
missile  are  parallel  to  each  flow  component 
direction  and  are  proportional  to  the 
directional  dynamic  pressure.  For  other 
shapes,  flow  field  similarity  in  the  crossflow 
regime  as  the  angle  of  attack  changes  is  the 
major  requirement  for  the  cross  flow  theory  to 
be  applicable.  Thus,  it  is  reasonable  to 
consider  extension  of  the  theory  to  sharp- 
edged  debris  missiles  that  force  boundary 
layer  separation  at  a  fixed  point  and, 
therefore,  produce  similar  potential  cross  flow 
fields  for  all  angles  of  attack.  In  principle, 
this  concept  allows  the  generation  of  lift, 
drag,  and  side  forces  for  certain  sharp-edged 
planar  symmetric  sections  if  the  drag 
coefficients  are  known  for  flow  normal  to  the 
three  major  faces  of  each  shape.  Normal  flow 
coefficients  can  be  found  in  the  literature  for  a 
variety  of  shapes.  The  final  form  of  the 
equations  includes  an  aspect  ratio  (tip  loss) 
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correction  for  finite  missile  dimensions,  and 
missile  face  porosity.  Table  2  illustrates  the 
form  of  the  equations  for  a  rectangular 
parallelepiped  shape.  Figure  2  compares  the 
cross-flow  predicted  lift  coefficients  for 
various  roll  angles  to  wind  tunnel  data  for  a 
parallelepiped  missile.  Summary  equations 
for  fourteen  basic  shapes  are  given  in 
Twisdale,  etal.  [1981]. 

RANDOM  TUMBLING  MODE  DRAG 

COEFFICIENTS 

Most  debris  shapes  are  not 
aerodynamically  stable  and  may  exhibit 
autorotation,  flat  rotation,  random  tumbling, 
coning,  or  other  motions  during  flight.  These 
motions  are  governed  by  the  initial 
conditions,  the  equations  of  motion,  and  the 
resulting  aerodynamic  and  non- aerodynamic 
forces  acting  on  the  missile.  Randomizing 
types  of  forces  may  arise  from  turbulence  and 
non-ideal  gas  flows,  explosive  products,  and 
missile  interactions,  such  as  debris-debris, 
debris-structure,  and  debris-ground  impacts. 
Hence,  random  tumbling  motion  drag 
trajectory  coefficients  are  often  used  in  2-D 
and  3-D  drag  models  to  predict  motion  of 
irregular,  bluff  body  shapes  and  fragments. 
Because  of  the  differences  in  the  random 
tumbling  mode  (RTM)  coefficients  presented 
in  the  literature  (e.g.,  Bates  and  Swanson 
[1967]  and  Redmann,  et  al.  [1976]),  a 
validated  equation  is  summarized  herein.  Use 
of  RTM  coefficients  in  2-D  or  3-D  drag 
models  should  be  used  only  to  get  the 
approximate  center  of  the  impact  dispersion 
pattern  and  not  to  estimate  the  debris 
dispersion  or  safety  distance. 

The  expression  for  the  expected  value  of 
the  drag  coefficient,  Cd,  of  a  tumbling  missile 
is 

r2n  r2  n  r»2n 

Cd=\  Cd (a, (5,  5) 

•J  0  */  0  J  0  //a  \ 

/  ( a ,  p,  S)  dddpda 


where  a,  P,  and  d  are  orientation  angles  as 

specified  in  Figure  1,  and  f  (a,  p,  S)  is  the 
joint  probability  density  function  describing 
orientation  likelihood.  For  a  cylinder  of 
diameter  d  and  length  L,  cross-flow  theory 
indicates  that 

Cd  (a,  p,S)  =  CD  (a) 

-  CDc  sin3  a  +  ^  CDa  Icos3  a I 

(3) 

Assuming  uniformly  random  spatial 
orientation,  f  (cc,  p,  5)  =  V8  tt2  sin  a  and  the 
expression  for  the  RTM  coefficient  is  derived 
as 

^=i(xc»r+fz:cD«)  <4> 

where  the  subscripts  a  and  c  refer  to  axial  and 
cross-flow  directions,  respectively.  This 
expression  yields  a  significantly  higher 
expected  value  than  the  previously  published 
results  of  Bates  and  Swanson  [1976]  and 
Redmann,  et  al.  [1976].  It  is  noted  that  this 
general  formulation  agrees  with  that  given  by 
Sentman  and  Niece  [1967].  Trajectories 
computed  using  the  random  orientation  model 
(drag  force  only)  with  high  update  frequencies 
are  shown  by  Twisdale,  et  al.  [1979]  to 
converge  exactly  to  the  impact  point  predicted 
by  2-D  trajectory  calculations  with  Cd  given 
by  Equation  4. 

SPIN-STABILIZED  TRAJECTORIES 

The  possibility  of  spinning  flat-plate  or 
disc-shaped  fragments  has  been  recognized  in 
the  explosive  safety  literature  (e.g.,  Moseley 
and  Whitney  [1980])  in  terms  of  probable 
maximum  debris  range.  High  in-plane  spin 
rates  (flat  rotation),  imparted  as  a  result  of  the 
explosion  effects  or  secondary  missiles  from 
failed  rotating  equipment,  can  lead  to  lifting 
forces  and  significant  out-of-plane  trajectory 
motion.  When  these  conditions  exist  in 
explosive  safety  problems,  they  should  be 
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TABLE  2.  CROSS  FLOW  AERODYNAMICSTFOR  PARALLELEPIPEDS. 


Missile  type 


Box,  beam,  plate,  frame 


Generic  shape 


Rectangular  parallelepiped 


Subcategories 

b/d  =  1 

l,d  *  3 

b/d  =  4 

b/d  -  10 

b/d  =  50 

Set  nos. 

5,6,18,19,20 

7,8,9 

10,11 ,21722,23 

12,13 

cDa 

2.05  Solid 
0.82  Frame/ 
truss 

2.05  Solid 

2.05  Solid 

2.05  Solid 

0.82  Frame/ 
truss 

2.05  Solid 

*-0b 

2.05  Solid 
0.82  Frame/ 
truss 

2.0  Solid 

2.0  Solid 

2.0  Solid 

0.80  Frame/ 
truss 

2.0  Solid 

CDd 

2.05  Solid 
0.82  Frame/ 
truss 

1.4  Solid 

1.0  Sol  i  d_ 

1.075  Solid 
0.43  Frame/ 
truss 

1 .575  Sol  id 

Skin  friction 
correction,  f 


1,  L/b  <  3 

0.41  +  0.59e-z(L/b-3),  3  <  L/b  <  4 
0.46  +  010061  (L/b) ,  L/b  *»4 


Aspect-ratio 

correction 


0.59  +  0r41e-20b/d 
0.59  +  0.41e-20dA 
0.59  +  0.41e-20bA 


CD 


(b/L)  Crja  ka  cos3o  +  Cnp  kp  |cos«  sina|p|  cos6|  sine.  + 
(b/d)  CDd  kd  si  n3  &  si  n3  a 


CL 


CS 


P 


-(b/L)  Coa  ka  |cosa|  cosa  sina  + 

Cob  kb  I c os 6  sinc.|p  |cos«l  cose  + 
(b/d)  Cob  kd  I  s i n3 6 1  sin2ocoso 

Cpb  kb  cos«|cos«|p-l  sin«  sinpa  - 
(b/d)  Cod  kd  sin6  |sin«|  cos6  sir/a 

0.5  +  0.15  (b/d)  +  1.35  (b/d)2 


Ref.  area,  A 


dL 


considered  since  they  will  influence  the  safety 
distance.  A  procedure  to  determine  if  the 
initial  in-plane  rotation  rate  is  sufficient  to 
stabilize  a  plate  or  disk-shaped  fragment  is 
summarized  herein  from  the  work  of 
Twisdale,  et  al.  [1984], 

For  purposes  of  developing  the  governing 
equations,  a  circular  disc  geometry  is 
assumed.  Coupling  between  the  resultant 
aerodynamic  force  vector  and  the  gyroscopic 
angular  momentum  vector  will  cause  the 
spinning  disc  to  slowly  precess  about  an  axis 
perpendicular  to  the  spin  axis.  The  minimal 


rotational  speed  required  to  maintain  this  spin 
stabilized  motion  and  the  resulting 
precessional  rate  is  developed  for  a  symmetric 
disk  rotating  about  its  center  of  gravity.  The 
disc  sector  rotates  about  an  axis  parallel  to  the 
Y-axis  and  passing  through  the  center  of 
gravity.  Referring  to  Figure  3,  let  the  Y-axis 
of  the  rotating  reference  system  lie  along  the 

axis  ofthe  disk.  Let  this  axis  make  an  angle  6 
with  a  fixed  vertical  axis  OB,  and  let  it 

precess^at  a  rate  Q  about  OB.  The  X-axis  is 
in  the  horizontal  plane  AOC  and  is 
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a.  Lift  Coefficients  for  8  =  0°,  60°  b.  Lift  Coefficients  for  8  =  30°,  90° 

Figure  2.  Comparison  of  Cross  Flow  Lift  Coefficients  to  Wind  Tunnel  Data  for 
Parallelepiped. 


Figure  3.  Conical  Precession  of  a  Disc  About 


b.  Aerodynamic  Center  of 
Pressure 

a  Vertical  Axis. 
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Hit 


perpendicular  to  both  OY  and  OB.  The  Z- 
axis,  perpendicular  to  OY  and  OX,  lies  in  the 
vertical  plane  through  OB  and  OY. 

Projection  of  Q.  on  the  three  rotating  axes  thus 
gives 


1953],  Equation  9  is  evaluated  from 
l’Hopital's  rule,  which  leads  to 


Nr 

2l y  0)y 


(10) 


Qx  =  0;  £2y  -  Q  cos  &,  Qz  =  Q  cos  6  (5) 

and,  since  OY  coincides  with  the  axis  of  the 
disk, 

ox  =  Qt=  0;  Oh  =  &z  -  sin  B  (6) 

The  angular  velocity  cay  of  the  disk  about  OY 
is  an  input  value  determined  by  the  initial 
conditions. 


Since  the  angle  0  will  always  be  less  than  or 
equal  to  nH,  so  that  cos  6  is  non-negative,  a 

steady  precession  (and  thus  a  constant  Q)  is 
possible  only  if  the  rotational  speed  of  the 
disk  isjsuch  that  the  radical  in  Equation  9  is 
non-negative,  i.e., 

($>j-[2NrIzcos6]m  (11) 


Euler's  dynamical  equations  of  motion 
reduce  to 

Mx  =  iJsin  6(Iy  C0y  -Iz  12  cos  0);  My  =  (T; 

Mz  =  0  (7) 

where  />  and  lz  are  the  moments  of  inertia. 
From  Equation  7,  the  maintenance  of  the 
assumed  motion  requires  only  a  moment 
about  the  X-axis.  This  moment  is  provided 
by  the  aerodynamic  force  normal  to  the 
circular face  of  the  disk  acting  at  the  center  of 
pressure,  CP  (Figure  3b).  Assuming  for  the 
circular  disk  that  the  center  of  pressure  is 
located  approximately  halfway  between  the 
"leading  edge"  of  the  disc  and  its  geomemic 
center,  the  moment  due  to  the  aerodynamic 
forces  is  estimated  by 

Mx  =  1/2  Nr  (8) 

where  N  is  the  aerodynamic  force  acting 
normal  to  the  plane  of  the  disc.  From 
Equations  7  and  8,  the  steady  precessional 
motion  maintained  by  this  moment  is 

^  l y  0)y±[ly  <a|  -  2 Nr  Iz  cos  0]1/2 
'"f-~  2/z  cos  e  (9) 

Neglecting  the  high  rate  of  precession 
corresponding  to  the  plus  sign  [Rauscher, 


Equation  1 1  suggests  that  coy*  may  be  as  low 

as  zero  when  6  =  90  degrees.  However, 
examination  of  Equation  10  shows  that  the 
rate  of  precession  corresponding  to  this 
minimum  rate  of  spin  is  infinite.  Thus,  while 

mathematically  there  is  no  minimum  (Oy 
required  to  spin  stabilize  the  disc  at  6=  K/2,  in 
practiceif  the  product  QlyOhy  is  not  sufficient 
to  absorb  the  moment  I/2A?r,  the  disc  will 
tumble  until  6  becomes  of  sufficient 

magnitude  that  coy  >  £Oy*,  at  which  point  the 
disc  will  become  spin  stabilized  and  will  then 
precess^  at  the  angular  velocity  given  by 
Equation  9. 

Table  3  summarizes  initial  angular 
velocities,  coy,  required  to  stabilize  several 
circular  disk  missiles  for  different  angles  of 
attack  (a),  initial  velocities  (v0),  and  inertial 

orientations  (6).  Three  steel  disks  were 
selected  with  weights  of  10,  100,  and  1000 
lbs,  and  radius/thickness  ratios  of  about  10. 

Solutions  are  shown  for  a  =  10  and  30 

degrees,  v0  =  100  and  500 ft! sec,  and  6  =  30, 
65,  and  90  degrees.  For  the  smaller  disk  (A), 
higher  C0y*  are  required  to  maintain  the  fixed 
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TABLE  3.  SPIN  RATES  REQUIRED  FOR  SPIN-STABILIZED  TRAJECTORIES  OF 
CIRCULAR  STEEL  DISKS. 


Initial  Conditions 

a  v0  6 

{deg)  (ft/sec)  (deg) 

0)y*  (rpm)  for  Disks  A,  B,  and  C 

ABC 

wt  =  10  lbs,  wt  =  100  lbs,  wt  =  1000  lbs, 

r  =  4.7  in  r  =  10.6  in  r  =  23.7  in 

thk  =  0.5  in  thk  =1.0  in  thk  =  2.0  in 

10 

100 

30 

54 

26 

12 

65 

38 

18 

8 

89 

8 

4 

2 

500 

30 

270 

128 

60 

65 

189 

89 

42 

89 

38 

18 

9 

30 

100 

30 

94 

44 

21 

65 

65 

31 

15 

89 

13 

6 

3 

500 

30 

468 

221 

105 

65 

327 

155 

73 

89 

66 

31 

15 

orientation  than  are  required  for  the  heavier 
disks.  For  all  three  disks,  minimum  spin  rates 

are  needed  as  6  — ►  90°  (horizontally  spinning 
disks).  The  required  spin  rates  increase 
significantly  as  v0  is  increased.  These 
example  calculations  indicate  that  for  low 
angles  of  attack,  and  horizontally  oriented 
disks,  less  than  one  revolution  per  second  (60 
rpm)  may  stabilize  >10  lb  disks  traveling  at 
less  than  500  ft/sec.  If  spin  stabilization  is 
possible,  trajectory  calculations  can  be  made 
using  the  RO  6-D  model  with  zero  update 
frequency  for  fixed  inertial  orientation  flight. 

MODEL  COMPARISONS  AND 

RESULTS 

The  random  orientation  trajectory  model 
has  been  developed  such  that  it  can  operate  in 
any  of  three  modes:  3-D  constant  drag, 
random  orientation  with  drag  force  only 


(random  drag),  or  random  orientation,  with 
drag,  lift,  and  side  forces  (full  random). 
Calculations  with  the  RO  6-D  model  as 
presented  for  the  three  secondary  debris 
missiles  summarized  in  Table  4.  The  ejection 
angles  in  Table  4  are  measured  from  the 
horizontal  axis. 

Steel  Joist.  In  ESB  Technical  Paper 
Number  13  [DoD,  1991],  the  steel  joist  in 
Example  Problem  1  is  estimated  to  travel  a 
maximum  distance  of  62  ft,  based  on  a 
MUDEMIMP  constant  drag  calculation. 
Using  the  RO  6-D  model  with  drag,  lift,  and 
side  force  components  and  an  update 
frequency  of  2  hz,  100  trajectory  simulations 
produce  the  impact  scattergram  shown  in 
Figure  4.  The  initial  position  of  the  missile  is 
X  =  0,  Y  =  0,  and  the  horizontal  component  of 
v0  points  in  the  positive  X  direction  The 
range  statistics  are  summarized  in  Table  5. 
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TABLE  4.  SECONDARY  MISSILE  DESCRIPTIONS  FOR  TRAJECTORY 
COMPARISONS. 


Parameters 


1.  Steel  Joist 


Missile  Description 
2.  Steel  Door 


3.  Secondary 
Fragment 


Weight  (lbs) 
d  (in) 

L/d 

V0  (ft /sec) 

EjectionAngle  (deg) 
Reference 


320 

8 

20 

111 

85° 

ESB  No.  13 
[DoD,  1991] 
Ex.  Problem  No.  1 


117 

36 

2.33 

677 

0° 

ESB  No.  13 
[DoD,  1991] 
Ex.  Problem  No.  2 


0.2 

0.53 

6 

985 

45° 

Kineke,  1976 


The  centroid  of  this  scattergram  is  indicated 
by  the  circle.  About  half  the  trajectories 
travel  further  than  the  mean,  and  about  20% 
exceed  the  62  ft  maximum  in  ESB  No.  13. 
Randonrtumbling  flight  with  drag  forces  only 
has  also  been  simulated  with  the  RO  6-D 
model  by  using  an  update  frequency  of  100 
hi.  The  predicted  impact  point  is  very  near 
the  centroid  of  the  data. 

Steel  Door.  The  results  for  the  steel  door 
l  DoD,  1991]  indicate  a  maximum  range  of 
808  ft.  -The  RO  6-D  scattergram  for  5  hi 
update  frequency  is  shown  in  Figure  5,  and 
the  range-to-first-impact  is  summarized  in 
Table  5l  The  RTM-predicted  range  is  226  ft 
to  first  impact.  These  results  are  sensitive  to 
update  frequency  with,  slower  updates 
yielding  larger  variances  in  both  x  and  y 
directions.  Higher  update  frequencies 


converge  to  RTM  results  with  minimal 
variance.  Note  the  cluster  of  points  at  about 
=  100  ft.  These  indicate  the  position  of  first 
impact  for  initial  orientations  that  result  in  a 
net  downward  aerodynamic  force,  which 
would  cause  a  skid  or  ricochet,  followed  by 
an  upward  rebound  of  the  door. 

Secondary  Fragment.  The  steel 
cylindrical  fragment  comparisons  are  given  in 
Figure j5,  and  the  statistics  are  summarized  in 
Table  5  for  an  update  frequency  of  10  hz. 

Similar  to  the  other  missiles,  the  fi  +  lcrand  fJ. 

+  2<rlclistances  significantly  exceed  the 
predicted  RTM  drag  range.  However,  the 
dispersion  pattern  is  significantly  smaller, 
reflecting  the  reduced  Lid  over  the  beam 
missile  and  the  higher  update  frequency  used 
for  this  smaller  fragment. 
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Figure  4.  Impact  Scattergram  for  Steel  Joist  Missile. 

TABLE  5.  RO  6-D  TRANSPORT  RANGE  STATISTICS. 

2-D  Drag  RO  6-D  Range  Distance1 


Steel  Joist 
Steel  Door 
Steel  Fragment 

t  RncpH  rm  1HO  cirr 


(M)  (o) 

52  15~ 

495  280 

1797  121 


fi+  la  ji  +  2a  Max 

67  82  ~  84 

775  1055  1083 

1918  2039  2252 
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X(ft) 

Figure  5.  Impact  Scattergram  for  Steel  Door  Missile. 


SUMMARY 

For  primary  fragments  and  secondary 
missiles  with  non-chunky  shapes  (plates, 
doors,  slender  fragments,  and  structural 
elements  with  high  slenderness  ratios), 
trajectory  models  that  treat  only  the  drag 
component  of  the  aerodynamic  force 
underpredict  maximum  debris  range  and 
dispersion?  The  RO  6-D  model,  which 
considers-  drag,  lift,  and  side  force 
components,  provides  an  efficient  alternative 
model  that  has  been  used  previously  in 
nuclear  power  plant  missile  risk  assessment 
and  to  recommend  debris  impact  velocities  for 
the  DOE  NPR  program.  An  aerodynamic 


library  for  typical  secondary  missile  shapes 
has  been  developed  for  use  with  the  RO  6-D 
model.  A  formula  for  evaluating  spin- 
stabilized  flight  potential  for  disk-shaped 
fragments  is  presented  and  evaluated  for 
several  disk  sizes  and  weights. 

The  trajectory  calculations  illustrate  some 
basic  features  of  the  RO  6-D  model  and  the 
fact  that  maximum  ranges  for  hazardous 
debris  density  may  be  underestimated  by 
current  prediction  methods.  The  model  would 
easily  be  validated  for  explosive  safety  siting 
analysis  through  a  series  of  calculations  with 
a  full  6-D  model  for  selected  shapes,  coupled 
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with  direct  statistical  analysis  and 
comparisons  to  Q/D  test  data.  The  prediction 
methodology  can  efficiently  simulate  both 
chunky  and  non-chunky  debris,  and  it  offers 
several  important  theoretical  advantages  over 
the  drag  models  for  non-chunky  shapes. 
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